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1 Introduction 

The Quadratic Assignment Problem (QAP) is the problem of making exclusive assignments of n 
indivisible entities to n other indivisible entities in such a way that a total quadratic interaction cost 
is minimized. The problem can be interpreted from a wide variety of perspectives. The perspective 
adopted in this paper is that of the generic facilities location/layout context, as in the seminal 
work of Koopmans and Beckmann [8]. Specifically, there are n facilities (or departments) to be 
located at n possible sites (or locations). The volume of traffic going from facility i to facility 
j is denoted fij. The travel distance from site r to site s is denoted drs- A quadratic "material 
handling" cost of hirjs is incurred if facilities i and j are assigned to sites r and s, respectively. In 
addition, there is a fixed cost (an "operating cost"), Oir, associated with operating facility i at site 
r. It is assumed (without loss of generality) that the units for "distance", "volume of traffic", and 
"operating cost" have been chosen so that the /ijrjs's and commensurable. The problem 

is that of finding a perfect matching of the facilities and sites so that the total material handling 
plus facilities operating costs is minimized. 

Let F := {1, 2, r/} and T := {1, 2, ...,<?} be the sets of facilities and sites, respectively. Without 
loss of generality, assume r] = q = n. For i £ F and r € T, let Wir be a 0/1 binary variable that 
indicates whether facility i is assigned to (or located at) site r {wir = 1), or not {wir = 0). Then, 
a classical formulation of the QAP is as follows (see [lOj ) : 
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Problem 1 (Problem QAP) 



minimize: 

(1) 

ieF jeF reT seT ieF reT 

subject to: 

^Wir= 1; reT (2) 



ieF 

Y,Wir= 1; ^Gi^ (3) 

reT 

u-irGjO,!}; ieF, reT (4) 

where . hi^-js — fij^rs ~1~ fji^sr (^) 

Problem QAP was shown to be NP-Hard as far back as the 1970's (see jl4|). Moreover, it 
has been known for some time that the Travehng Salesman Problem (TSP; see [S]) and other NP- 
Complete combinatorial optimization problems (see [7] , [11] , or [12] ) can be modeled as special cases 
of the problem. Hence, the thrust of research on the problem has been towards the development 
of heuristic procedures and "tight" lower bounds (see [I], [3], [10], and [13] for extensive reviews). 

In this paper, we present a new linear programming (LP) formulation of the Quadratic As- 
signment Problem (QAP). The modeling approach is similar to those in [5] and [6]. However, the 
proposed model is an order of size smaller, having 0{n^) variables and 0{n'^) constraints, where n 
is the number of assignments. A small experimentation that was undertaken in order to gain some 
rough indications about the computational performance of the model is discussed. 

The plan of this paper is as follows. The proposed linear programming formulation is developed 
in section 2. Computational testing and results are discussed in section 3. Conclusions are discussed 
in section 4. 



2 Development of the Formulation 

Our approach consists of modeling the QAP in the framework of the multipartite graph G = (V, 
A) illustrated in FigureUl Nodes in this graph represent (facility, site) pairs. Arcs are labeled with 
triples (a,6, c) S {F,T,F), and represent assignment decisions. 

Remark 2 The association we make between Graph G and the QAP is to interpret a positive flow 
into/out of any node of the graph to mean that the corresponding facility and site pair have been 
assigned to each other. 

Definition 3 

1. The set of all the nodes of Graph G that have a given facility index in common is referred to 
as a level of the graph. The i*^ level [i e F) is denoted Li{G) := {{u,v) eV\u = i}. 

2. The set of all the nodes of Graph G that have a given site index in common is referred to as 
a stage of the graph. The j*^ stage (j e T) is denoted Sj{G) := {{u,v) eV \ v = j}. 



2 



3. A path in Graph G that simultaneously spans the levels and the stages of Graph G is referred 
to as a perfect bipartite matching (p.b.m.) path of the graph. 

4. The set of all the p.b.m. paths of Graph G is denoted CI. That is, 

n := {{{ii, 1, i2), {i2,2, is), {in-i , n - 1, in)) e ^"-1 : V 7^ ^g V {p, q) G (T, r\{p})}. 



Figure 1: Illustration of Graph G 

Remark 4 

1. There exists a one-to-one correspondence between perfect matchings of the facilities and sites 
and p.b.m. paths of Graph G; 

2. There exists a one-to-one correspondence between the perfect matchings of the facilities and 
sites and feasible solutions to Problem QAP; 

3. There exists a one-to-one correspondence between the p.b.m paths of Graph G and the 
feasible solutions of Problem QAP. 

Assumption 5 Throughout the rest of this paper, it is assumed (w.l.o.g.) that: 

1. The number of assignments is greater than 5 (i.e., n > 5); 

2. All vectors of variables are column vectors. 

Notation 6 The following notation will be used throughout the rest of this paper: 
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1. The set of real numbers is denoted by TZ; 

2. For two column vectors a and b, ^ ^ ^ = (a^,b^)^ will be written as "(a, b)" (where 
(•)^denotes the transpose of (•)), except for where that causes ambiguity; 

3. The z*'* component of a column vector a is denoted Eij; 

4. The notation "0" denotes a column vector of comfortable size that has every entry equal to 
0; 

5. The notation "1" denotes a column vector of comfortable size that has every entry equal to 

1; 

6. The convex hull of (•) is denoted Conv{-); 

7. The set of extreme points of (•) is denoted Ext{-); 

8. The set of stages of Graph G from which arcs of Graph G originate is denoted R\ i.e., 
R:=T\{n}; 

9. V {i, j,u,v, k,t) € F^, M (p,s) & B? : 1 < p < s, Zi^i^jupvkst denotes a non-negative variable 
that represents the amount of flow in Graph G that propagates from arc (z,l,i) on, to arc 
{k,s,t), via arc {u,p,v); 

10. V ii,j,k,t) G F^, V (r,s) G i?2 . 1 < ^ <- 

Uirjkst denotes a non-negative variable that 
represents the total amount of flow in Graph G that propagates from arc {i,r,j) on, to arc 
{k,s,t). 

2.1 Model Constraints 

The constraints of the Integer Programming (IP) version of our proposed model are as follows: 

^iA,n<^,vv,3,t = 1 (6) 

ieF jeFveFteF 

vGF vGF 

i,j,k,t,u€F;p,s€R:l<s<p<n — l (7) 

veF veF 

i, j, k,t,u G F] p, s £ R : 1 < p < s — 1, s > 3 (8) 

Vijljupv ^ ^ ^ ^ Zi,l,jupvkst ~ 0) 
kGF t£F 

i,j,u,v€F]p,seR:l<p<s (9) 
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yi,l,jkst ^ ^ ^ ^ ^i,l,jupvkst — 0, 

i,j,k,tGF; p,seR:l<p<s (10) 



p€R:l<p<s veF peR:p>sv€F 

i,j,k,t,ueF; seR\{l} (11) 



Uupvkst ^ ^ ^ ^ ^i,l,jupvkst 0; 
i£F jeF 

u,v,k,t£F; p,sGR, l<p<s (12) 



E yirjkrt + E yirjk,r+l,t = 0] iJ^F; r£R (13) 

ik,t)£F2 (k,t)eF'^: 

k=/=j; {k,r+l,t)eA 



E EEy--+ E EE«n E E^^n- + 

seR: s>r keF teF sGR: s<r keF teF seR:s>r+l kGF 

~l~ ^ ^ ^ ^ Uirjisk + ^ ^ ^ ^ Uirjksj ~l" ^ ^ ^ ^ Uirjjsk ~ 0; 
sGR:s>r+l fceF seR:s>r+l fceF seR:s>r+2 keF 

i,jeF; reR (14) 
yirjfcst G {0, 1}; i,j,k,teFr,seR:r<s; (15) 
Zi,i,jupvkst G {0, 1}; i, j,u, V, k, t e F; p,seR: l<p<s (16) 



Constraint (l6|) initiates the propagation of one unit of flow from stage 1 to stage 3 of Graph G. 
Constraints ([7]) stipulate that the total amount of flow from arc (i, that propagates through arc 
{k,s,t) and enters node {u,p + 1) is equal to the amount of flow from arc (z, l,j) that propagates 
through arc {k,s,t) and leaves node {u,p + 1). Constraints ([8]) stipulate that the total amount of 
flow from arc (i, 1, j) that enters node {u,p+l) to propagate on to arc {k, s, t) is equal to the amount 
of flow from arc (z, l,j) that leaves node {u,p + 1) to propagate on to arc {k,s,t). Constraints ([9]) 
and (jlOp ensure that the propagation of the flow from a given arc at stage 1 of Graph G onto a 
given arc at another given stage of the graph is consistently accounted across all the other stages 
of the graph. Constraints ()12p stipulate that the total amount of flow that propagates from arc 
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{u,p,v) onto arc {k,s,t) is equal to the total of the flows from arcs at stage 1 that propagate onto 
arc {k,s,t) via arc {u,p,v). Constraints (|lip stipulate, essentially, that the total flow on any given 
arc of Graph G must propagate on to every level of the graph, or be part of a flow propagation 
that spans the levels of the graph. Constraints (jl3p ensure that the initial flow propagation from 
any given arc occurs in an "unbroken" fashion. Finally, constraints (|14|) stipulate (in light of the 
other constraints) that flow from arc {i,r,j) of Graph G cannot propagate back onto neither level 
i nor level j of the graph. 



Theorem 7 The number of variables in the system is 0{rfi); The number of constraints 

in the system ([^-|77P is 0{n'^). 

Proof. Trivial. ■ 



Definition 8 



1. We refer to the set of points in the space of the y- and z- variables that satisfy the system 
(f6l)- (fT6]) as the "IP Polytope," and denote it by Q/; i.e., Qi := {{y,z) € TZ^ : (y, z) satisfies 
([6|)- ()16p }. where ^ is the number of variables in the system ([6l)- (|16p . 

2. We refer to the linear programming relaxation of Qj as the "LP Polytope," and denote it by 
Ql\ i.e., Ql := {{y,z) € TZ^ : {y,z) satisfies (j6|)-p^. and < (y, z) < 1}, where ^ is the 
number of variables in the system (|6l)- (ll4|) . 

3. The set of points in the space of the tt;-variables that satisfy the system (l2|)-(|4]) is referred 
to as the "Assignment Polytope," and is denoted by Wj; i.e., Wj := {w £ TZ'^ : w satisfies 
©-(II])}, where u is the number of variables in the system (l2])-([4l); 

4. The linear programming relaxation of Wj is denoted by Wl] i.e., Wl := {w £ TZ^ : w 
satisfies ([I])-©, and Q <w < 1}, where v is the number of variables in the system ©-([Ij). 

2.1.1 Structure of the IP Polytope 

Theorem 9 (y, -z) € Qi <^=^ 3 exactly one set of facility indices, {v £ F, r = l,...,n}, such 
that: 



Uarbcsd 



and 



^a,l,bcrdesf 



1 for r, s £ R : r < s, and 
{a,b,c,d) = {ir ,ir+i,is, is+i ] 

otherwise 



1 for r, s £ R : I < r < s, and 

(a, b, c, d, e, f) = {ii,i2,ir,ir+i,'is,'is+i] 

otherwise 



(17) 



(18) 
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Proof, a) Let {y,z) E Qj. Then, given ([T5]), and (USD: 
Constraint ^ =^ 3 a unique set of facility indices, {ir € F, r = 1, . . . , 4}, such that: 

^il,l,i2i2. 2,1343, 3,14 ~ 1 (19) 

Condition (jlSp follows directly from the combination of (|19p. ([Tj), and ([8]). 

Condition (jl7p follows from the combination of Condition (jlSp with constraints ([9])- (jl0p and ()12p . 
6) <;^: Trivial. ■ 

Theorem 10 There exists a one-to-one correspondence between the feasible solutions to the system 
(^-(I^, and the perfect matchings of the facilities and sites. 

Proof. Combining Theorem^ with, constraints (jl4p . (jlSp . and (|16p . we must have: 

(i) p!7P and ([TH]) are satisfied for {y,z), and 
(20) 
ii) ip/ig V(p,g) G (T,r\M) 

The theorem follows directly from the combination of ()20|) , the definition of the y- and z- variables 
(see Notation [6l9 and[6ll0). Definition [3l 4, and Remark[Hl. ■ 

Corollary 11 Qj is isomorphic to 0, and to Wj, respectively. 

Definition 12 We refer to the perfect matching of the facilities and sites correponding to {y, z) G 
Qi as the "assignment corresponding to {y,z)," and denote it by the ordered set M{y,z) := 
(ii,i2, • • • ,in)) where iq is the index of the facility assigned to site q in the matching. 

2.1.2 Structure of the LP Polytope 

Lemma 13 (Flow propagation lemma 1) Let {y,z) E Ql- The following holds true: 

V(il, i2, ^3) ^4) € F^, 2/21,1,12,43,3,14 > <;=^ Zi^,l,i2,i2,2,i3,Vi,3M > 0- (21) 

Proof. Using constraints ^ and (I13p . constraints (llOp for p = 2 and s = 3 can be written as: 

?/ii,l,i243,3,j4 — ^ii, 1,4242, 2,4343,3,44 = V(ii, ^2, ^3) ^4) € F^ (22) 

The lemma follows directly from (j22p . ■ 

Lemma 14 (Flow propagation lemma 2) Let {y,z) G Q^. Then, we must have that: 
yr £ R:r>A, V(ii, i2, is, Wi) G F^, 

i) ^41,1,42,42,2,43,43,3,44 ^ 0;' 
^41, 1,42,43, 3, 44,4r,r,4,+l ^ ^ ^ ^ (^3) 

ii) '^4l,l,42,42,2,43,4r,r,V + l ^ 0. 
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Proof. 

a) Condition i. Using constraints Q, 

Zh,i,i2,i3,3,i4,ir,r,ir+i > =^ i , 1 ,13 ,3,24 > Vr G i? I T > 4, V(ii , 42 , «3 , ^4 , V , V+1 ) e -^'^ (24) 
From Lemma [131 

yh,i,i2,i3,3M > 2;^^, 1,^2,^2,2,43,43,3,24 > y{ii,i2,i3,u) ^ F'^ (25) 
Condition i) follows directly from (j25p . 

6) Condition ii. Using ([9|), (jlOp . (|12p . and (jl3p . constraints ([8]) for p = 2 and u = i^, can be written 
as: 

2ii,i,j2,j2,2,j3,is,s,j^+i - ^ 2;j-^,i,j2,i3,3,2;,js,s,i,+i = V(n, i2, ^3, ^s, ^s+i) ^ Vs E : s > 4 (26) 

Hence, in particular, we must have: 

Zii,i,i2,i2,2,i3,ir,r,ir+i — ^11,1,12,43, 3,M,ir,r-,v+i ^ y{ii,i2,i3,U,ir,ir+i) G V r G i? : r > 4 (27) 
Condition ii follows directly from (j27p . ■ 

Lemma 15 (Flow propagation lemma 3) The following holds true for all {y,z) Ql: 
\/r e R : 2 < r < n - 3, V(ii, i2, V, V+i, V+3) ^ F^, 

(i) Zij^,l,i2,ir+l,r + l,ir + 2,ir+2,r+2,ir + 3 > Oj 
(28) 
^^) ^«l,l,i2,ir,r,jr + l,ir+l,r' + l,4r+2 > 0- 

Proof, a) Condition i. Using ([9|), ([10]), ([T2|) . and p^ . constraints ([8]) for p = r, s = r + 2 , 

and n = ir+i, can be written as: 

^il,l,i2,f,r-,ir.+ l,ir + 2,'' + 2,V+3 ~ ^4l,l,42,4r+l,»' + l,ir+2,ir + 2,T' + 2,ir-+3 ~ 

VrGi?:2<r<n-3, V(ii, ^2, V+i, V+2, V+3) G (29) 
Hence, in particular, we must have: 

^il,l,i2,jr,»',4r+l,*r+2,f + 2,V+3 ~ , 1 ,«2 , V+1 ,»'+l,«r + 2 , V+2 ,''+2,ir + 3 — 

Vr€i?:2<r<n — 3, V(ii, ^2; V+ij V+2) V+3) ^ -^^ (30) 
Condition i) follows directly from (j30p . 

6) Condition ii. Using constraints ([9]), ([TO]) . ([T2]) . and ([T3]) . constraints ([7])forp = r + l,s = r, 
and « = ir+2, can be written as: 

Zii,l,i2,ir,r,ir+l,ir+l,r + l,ir+2 ~ ^ ^ ^il,l,i2,jr,n*r + l,4r+2,r' + 2,D = 

VrGi?:2<r<n-3, V(ii, i2, ir+i, V+3) £ (31) 



Hence, in particular, we must have: 

Zii,l,i2,ir,r,ir+l,ir + l,r + l,ir+2 ~ ^ii,l,i2,ir,r,ir+i,ir + 2,r + 2,ir+3 — 

\/r e R : 2 < r < n - 3, V(ii, 22, ir, V+i, «r+2, V+s) G (32) 
Condition ii) follows directly from (j32]) . ■ 

Notation 16 For {y,z) G Ql: 

1. The sub-graph of G induced by the positive components of {y,z) is denoted as: 

H{y,z): = {P{y,z)),E{y,z)), (33) 

where: 



(i,l) GF: ^^yi,ijj,2,t > Oiu 
jeFteF J 

{i,r) eV J2 Va^iMrj + 

a&F beF jeF:{i,r,j)&A 

+ Z^IZ ya,l,bj,r-l,i > 

aeF beF jeF:{j,r-l,i)£A 



\ 



(34) 



Eiy,z) ■.= ai,l,j) e A:^yi^ijj^2,t > ol 
I teF ) 



U 



(i,r,j)G^:^^ ^ ya,i,birj > 0} ■ (35) 

2. The set of arcs of H{y, z) originating at stage r of H{y, z) is denoted Tr{y, z); 

3. The number of arcs originating at stage r of Graph H{y,z) is denoted ')j.{y,z) = \Tr{y.,z) \ . 
For simplicity jj.{y,z) will be henceforth written as 7^ (unless that causes ambiguity); 

4. The index set associated with Tr{y,z) is denoted Ar{y,z) := {1, 2, . . . , 7^}. For simplicity 
Ar(y, z) will be henceforth written as A,.; 

5. The arc in rr{y,z) is denoted as ar^^{y,z). For simplicity ar,^{y,z) will be henceforth 
written as ar^u', 

6. The tail of Ur^u is labeled ir^uiy,z); the head of ar^y{y,z) is labeled jr^u{y,z). For simplicity, 
ir,uiy, z) will be henceforth written as v,^) and jr,u{y, z), as jV,^; 
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7. Where that causes no confusion (and where that is convenient), for r G R with r > 2, and 
(a, p) G (Ai, Ar) "yn,c,,i,ii,c,v,p,r,j/V,p" wih be henceforth written as "y(i,a)(r,p)-" Similarly, for 
(r,s) G with 1 < r < s and (a,p,cj) G (Ai,A^,A5), "zii,„,i,ji_„,i,_p,r.jv,p,j,,„,sj,,<," wih be 
henceforth written as "2(i,a)(r,p)(s,o-)-" 

Lemma 17 (Flow conservation lemma 1) Let {y,z) G Ql. The following holds true: Va G Ai, 
V(p, g, r, s) G i?^ : 1 < p < 1 < p < r < s, 

E E 

Proof. Va G Ai, y{p, q,r, s) ^ : 1 < p < q; 1 < p < r < s, 
Yl Yl ^(i,»)(P,upK<],uq) = Yl y(i-"){p-^p) (Using ©) 

VpGAp UqGAq UpeAp 

= Y Y H'^,cKp,'^p)(r,i^r) (Using dSD) 

UpeAp UreAr 

= Y y{ha){r,ur) (Using ([10])) 

= Y Y H'^,c.)(r,ur)is,u,) (Using Q) 

UreAr UseAs 

m 

Definition 18 {''Paths in (y,z)''^) Let{y,z) G Ql- V (r, s) G i?^ : s > max{3, r+1}, V(z^i, z/^, z^s) G 
(Ai,Ar,As), a set of arcs of H {y , z) , 

{iar,ur,---,as,uj G (-E(y,z))'^-''+i : > V(p, g) G i?^ : 

max{2, r}<p<g<s-l; ip^i.^ = jp-i,,yp_-^ Vp G (i? n [r + 1, s])} 

is referred to as a "path in {y,z) from (r, z^^) to (5,1/^)." 

Notation 19 Let {y,z) G Ql- V (r, s) G i?^ : s > max{3, r + 1}, V {p,a) G (A^, A^), 

1. The set of all paths in {y,z) from {r,p) to {s,a) is denoted U(r,p)(s,cT){yT z)] 

2. The index set associated with C^(r,p)(s,o-)(y) -2) is denoted ^i^r,p){s,a){yj := 
{1,2,..., '^(^r,p){s,a) (y, z)}, where V'(r,p)(.,a) (y, ^) := I U(r,p)(s,cT) {y, z) I ; 

3. The k^^ element of C/(r-,p)(s,<7)(y, ^) ^ ^(.r,p)(s,a)(y, ^)) is denoted C^j.^p^^(s,a),k{y^^)- 

Theorem 20 (Distinct paths in (y,z)) Let {y,z) G Ql- Then, \/ {r, s) ^ R'^ : s > max{3,r+l}, 

V(ai,a2) G A2, V(/3i,/32) G A^ : ?7(^,a^)(^^^^)(y, z) 7^ ^■,U^,.,a2){s,f32)iy^ z) / 0, V /c G ^(r,ai){s,/3i)(y, 2=), 

V t G ^>(r,a2)(s,/32)(y, 'C(r,ai),(s,/3i),fc(y, 2) 7^ '^(^^^2)^(5,^2).* ^) 3 {5 G i? : r < 5 < s; (7i, 72) 

G(A„A,\{7j)}9{ e ^(r,a2),(s,/32),t(y>^)} • 



Proof. The theorem follows directly from the combination of constraints ([9]), (jlOl) . ()12p . and 
(fT3|) . and Definition [HI ■ 
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Theorem 21 (Path structure theorem 1) Let {y,z) E Ql. The following holds true: 
V(r,s) G with s > max{3,r + 1}, V(p, cj) G (Ar,A^), y(r,p){s,a) > <^ U(^r,p){s,a)iy, z) + 0- 

Proof. First, note that it follows directly from the combination of Lemmas 1131 [T^ andll5|, that 
the theorem holds true for all (r, s) G i?^ with s G {r + 1, r + 2}, and all (i/^) ^s) S (Ar, Ag). 

a) ^=^: Assume there exists an integer w > 2 such that the theorem holds true for all (p, t) 
G with t = p + uj, and all (z^p, z^f) G (Ap, At). We will show that the theorem must hold for all 
(p, u) G B? with u = t + 1 = p + uj + 1, and all (z^p, Vu) G (A^, A„). 



Let (p, u) G i?^ with n = p + a; + 1, and (i^p, i/^) S (Ap, A^) be such that: 



> 0. 



Define: 



B(p,up)iu,u^){y,z) ■■= {a e Ai : Z(i^a)ip,up)(u,u^) > 0}. 
Then, P). and ([36])^ 

i) B(p ){y,z)^0, with 



(36) 
(37) 

(38) 



Condition (j38p . and constraints ([8|) and ()13p 

V a G B(p,j.p)(„,^„)(y,2;), 3 C„,(p,^p)(„,^^)(i/, z) C Ap+i 9 : 

i) ip+1,/3 = jp,,.p V ^ G Ca,(p,i.p)(„,,.„)(y,z) 

){y,z); and 

Z{l,a){p,Up){u,Uu) - Yl ^(l,a){p+l,/3)(«,i/„)- 

By assumption (since u = {p + 1) + w), condition ({391 ^) =^ 

{y,z)^% yi5eC^ 

,{p,Up){u,Un) iy,z) (40) 
Also, it follows from the combination of condition (j36p . constraints ([7]), and constraints (jlip . that: 

V /? G C„,(p,^p)(„,^^)(y,z), 3 {T(p_^p)(p+i^^)(„_^^)(y,z) C $(p+i^^)(„_^^)(y,2;)} 3 : 
{(^)(i,Q)(p,i.p){g,i.,,0 > G T(p^j,^)(p+i^^)(„^^^)(y,z), 

Vg G (i?n [p+ l,it]), and V i/,,, G A, : Op^j.^^^ G /^(p+i,/3)(n,i.„),i(y, 2)} (41) 
Hence, V /? G C„^(p_^p)(„^j,^)(y, z) and V t G T(p^^,p)(p+i^^)(„_^^)(y, z), 

L ■■= {^{p+i,l3){u,u^)Ay^^) U {ap,i.p}) 
is a path in (y,z) from {p,i^p) to {u,Vu)- Hence, we have that C^(p ^.^^(^ j^^)(y, z) 7^ 0. 
b) <^^: Follows directly from Definition psp and constraints p2p . ■ 



(39) 
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Corollary 22 Let {y,z) E Ql. The following hold true: 

i)y S e R\{1} , V(a,cr) G {Ai,As), y{l,a){s,a) > 7^ 

V (r, s) e (ii\{l})2 with s > max{3,r + 1}, V(a,/9,cj) G (Ai, A^, A^), 

ii.l) U(^i,a){s,a)(.y,z) + 0, and 



^(l,Q){r,p)(s,(T) > 



Definition 23 ("p.b.m. path in (y,z)") Let {y,z) G Ql- V(z^i,i^„_i) G (Ai,A„_i), a path in 
from (1, z^i) to (n — l,zv„_i) is referred to as a "perfect bipartite matching (p.b.m.) path in 
{y,z) (from (1, i^i) to (n,z^„_i))." 

Notation 24 Let {y,z) G Ql- For all (a,/3) G (Ai,A„_i), 

1. The set of ah paths in {y,z) from to [n — 1,(3) is denoted as Yiap{y^z); 

2. The index set associated with Ilap{y,z) is denoted '^af3[y-,z) := {1, 2, ... , n ^piy , z)} ., where 

T^apiy^z) := \Uap{y,z)\; 

3. The /c*^ element of Ilaisiy, z) is denoted Va^kiy, z). 
Remark 25 Let {y,z) G Ql- V {a, 13) G (Ai,A„_i), 

2. 2) („_l,/3)(2/,z); 

4. We assume (w.l.o.g.) that: Vai3k{y,z) = J^{i,a),{n-i,f3),k{y^ ^) ^ ^ £ ip^p{y,z). 

Theorem 26 (Equivalence of p.b.m. paths and 2-matching solutions) For{y,z) G Ql, ev- 
ery p.b.m. path in {y, z) corresponds to exactly one perfect matching of the facilities and sites. 

Proof. Definition 1181 constraints ()14p . and Defintions [3l3 -[3j4 imply that every p.b.m path in 
{y,z) corresponds to exactly one p.b.m path of Graph G. The theorem follows directly from the 
combination of this Remark [HS. ■ 

Theorem 27 ("Convex independence" of p.b.m. paths) For {y,z) G Ql, a given p.b.m. 
path in {y,z) cannot be represented as a convex combination of other p.b.m. paths in {y,z). 

Proof. Theorem 1261 implies that every p.b.m. path in (y, z) corresponds to an extreme point of 
Conv(WL) = Conv{Wi). The theorem follows directly from this. ■ 

Theorem 28 (Path structure theorem 2) Let {y,z) G Ql- The following holds true: \/ r £ R, 
V p G Ar, 3 {(a,/3) G (Ai, A„_i); t G ^'Q,/3(y,z)} 3 ar,p e Va/3iiy, z). 
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Proof. 

Case 1: r = 1. From ([Ml) and ([35]) : 

ar,p e E{y, z) ^ 3 a G A2 3 y(r,p)(2,a) > 0. 
Condition (j42p and constraints ([9]) 

3 7 G A„_i 3 Z(^r,p){2,a){n-l,'y) > 0. 

Condition (j43p and constraints ([9|) 

3 7 G A„_i 3 y(r,p)(n-l,7) > 0. 

The theorem follows from the combination of ()44p with Theorem 1211 
Case 2: r = n-1. From ([HI) and ([35]) : 

ar,p e z) ^3 a G Ai 3 y(i,a){r,p) > 0. 
The theorem follows from the combination of ()45p with Theorem 12 1[ 
Case 3: 1 < r < n-1. From pi) and 

ar,p e z) ^ 3 a G Ai 3 y(i,a){r,p) > 0. 
Condition (j46p and constraints ([9]) 

37 G A„_i 3 2(l,a)(r,p){n-l,7) > 0. 

The theorem follows from the combination of (j47p with Corollary I22^i. ■ 

Corollary 29 Let {y,z) G Ql- The following hold true: 
i) Va G Ai, Vr G ii\{l}, Vp G A^, 

2/{i,a)(r,p) > 3 (/3 G A„_i; i G z)) 3 a^,p G VaPiiy, z); 

a) V(r, s) G -R^ : 1 < r < s, V(a, p, a) G (Ai, A^, A^), 

^(i,«){r,p)(s,a) > <^ 3 (/3 G An-i; i G ^0/3(^,2)) 3 (o^.p, fls.a) G 7^^^,(y,2;). 
Lemma 30 (Flow conservation lemma 2) Xei (y, z) G Ql. T/ie following hold true: 

i/pGAp /3eA„_i t6*a/3(?/,z): 

Va G Ai, V(p, r) £ : I < p < r, \/ Vr ^ K 
Va G Ai, V(r, g) G i?^ : 1 < r < g, M Ur ^ K 
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Proof. The lemma follows directly from the combination of constraints ([9]) and (jlOp . and 
Theorems [201 [271 and[28l and Corollary [ 



Definition 31 (p.b.m. path "weight") Let {y,z) E Ql- For (a, (3) G (Ai,A„_i) such that 
y(i,o)(n-i,/3) > 0, and k G z), we refer to the quantity 

(p,g)G-R^; {i^p,i^q) e (Ap.Aq): 

l<p<g; (ap,,.p, a,,^,) G P^i^^,{y,z) 

as the "weight" of (p.b.m. path) Vapk{y^z). 

Remark 32 It follows directly from Definitions and\3^ that for {y,z) G Ql, ^ai3i{y,z) > V 
(a,/?) G (Ai,A„_i) : ^^piy^z) / 0, V 6 G Po,(3,{y,z). 

Theorem 33 (Path structure theorem 3) Let {y,z) G Ql- The following hold true: 
V a G Ai, V r G R\{1}, V p G A^, 

2/{l,a)(r,p) = Yl Yl UJal3c{y,z) 

/3eA„-i tg<I'^^(j/,2): 

a) V(r, s) G i?^ : 1 < r < s, V(a, p, a) G (Ai, A^, Ag), 

^(l,a)(r,p)(s,cr) = ^af3Liy,z) 

Hi) V(r, s) G i?^ : 1 < r < s, V(p, a) G (A^, A^), 

2/{r,p){.,<7) = X] Yl ^o^Piiv^z) 

aGAi/3GA„_i tG'I'a,3(j/,2:): 

Proof, a) Condition i. First, note that from the combination of constraints ([6]), ([7]), ([8]), and 
(jlSh : Remark I32| and Theorems (I27p . and (I28p . we must have: 

Y Y Y Hl,c.)i2,m,5) = Y Y Y ^aAy^z) = l. (51) 

oGAi /3gA2: <5gA3: oGAi eGA„_i i.e'^ag{y,z) 

i2,l3=jl,a h,5=j2,f3 

a.l) From Definition [311 dSIl) ^ 

X] Y ^(l,")(2,/3)(3,5) = X] ^ocgiiy^z) Va G Ai (52) 

/3GA2: (5gA3: eGA„_i Le^ae{y,z) 

Lemma [T7] and relations (j52p 



X] X] " X] X] '^aet(y'^) Va G Ai,V(r,s) G i?^ : 1 < r < s (53) 

pGAr ctGAs fSAn-l tG'I'QeCj/.z) 
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Using constraints 1^, ()53p can be re-written as: 



y{l,a)ir,p) = ^agc{y,z) Vo G Ai, Vr € -R\{1, 71 - 1} 

Using Theorem [20l (j54p can be written as: 

^ y(i,a)(r,p) = 5Z ^c^Av^z) Vq E Ai,Vr G i2\{l,n - 1} 

Or.p 6 'PagL{y,z) 



Re- arranging (j55p gives: 



E 



\ 



V 



Va € Ai,Vr e i?\{l,n 



a. 2) Combining Lemma [301m with Definition 1311 we have that: 



SEE E 



ar,p G VagL{y,z) 

V r e ii\{l, n - 1}, yqeR:q>r, V(a, /9) S (Ai, A^ 
Relations (1561) and (IFT 



2/{i,a)(r,p) = X] X] Vr G i?\{l , n - 1} , V(a, p) G (Ai , A^) 

eGA„_i tGl'aeCy.z): 

a. 3) Using constraints (jlOp . (j53p can be re- written as: 

X] y(i.")(^><^) = Y J2 ^a(,c{y,z) Va G Ai,Vs G i? : s > 2 

o-GAs eGA^-i ;,G*ap{j/,2) 

Using Theorem [201 dSSD =^ 

X] = X] X] X] '^aei(yi^) Va G Ai,Vs G : s > 2 

o-GAs (tGAs eGA„_i tG'I'ae{j/,2;): 
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Re- arranging (j60|) gives: 

/ 



E 



\ 



as,cTe'Pa£it(3/,2;) 



Va G Ai,Vs € i? : s 



a. 4) Combining Lemma [30li with Definition 13 H we have that: 

2/(l,a){s,<x) = X] Yl Hho^)ip,Up){s,a) 

2 E E E 

V(p, s) e -R^ : 1 < p < s; s > 2, V(a, cr) G (Ai, A^) 
a. 5) Relations (lUTD and (P]) ^ 

= X] X] ^a£,t(y,z) Va G Ai,Vs G : s > 2 

egA„_i (,evi'Qp(j/,2): 

0.6) Condition i of the theorem fohows from the combination of ()58p and (j63p . 
6) Condition ii. 

6.1) Using Theorem 1201 and Corollary [291m, ()53p can be re-written as: 

peAr o-eAs /3eA„_i ie>i'„/3(y,z): peAro-eAs 

(ar,p, as,<T)eP2^^(?/,2) 

V(r, s) G i?^ : 1 < r < s, Va G Ai 
Re-arranging ([M|) gives: 

f \ 



peAr creAs 



Z{l,a){r,p){s,a) - Z Z ^0(3^11, z) 



V(r, s) G -R^ : 1 < r < s, Va G Ai 
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6.2) From Lemma [30ln, we have: 

y{l,a){s,a) = ^(l,a)(r,p)(s,a) + X] X] X] ^{l,a)(r,i/,)(s,<7) 

V(r,s) G i?2 . ;l < r < s, V(a,p,cj) G (Ai,Ar,A^). (66) 

6.3) From Condition i, we have: 

(a,,p, as,c.)GP2^,(s/,^) '^'-^^ (a,,,,,a,,<,)GP2^,(y,2) 

V(r,s) € -.Kr < s, V(a,p,cj) G (Ai,A^,A,) (67) 

6.4) Definition EI] 

V(r, s) G i?^ : 1< r < s, V(a, p, a) G (Ai, A^, A,) (68) 

6.5) Relations (IMD - 



(ar,p,a3,a)eP^^^(?;,2) 

V(r, s) G i?^ : 1< r < s, V(a, p, a) G (Ai, A^, A,) (69) 
6.6) Combining (|65p and (|69p . we must have: 

2;(l,a){r,p)(s,z.,) = X] X] ^aPiiy^z) 

V(r, s) G i?^ : 1< r < s, V(a, /J, cr) G (Ai, A^, A,) (70) 
c) Condition Hi. From the combination of constraints (1121) and Condition ii, we have: 



y{r,p){s,a) = Y = Yl Y Y UJaMV^z) 

(ar,p, as,<T)eP2^^fe,2) 

V(r, s) G i?^ : 1< r < s, V(p, a) G (A^, A,). (71) 



Corollary 34 (y, z) £ Ql <^=^ (2/)-^) corresponds to a convex combination of perfect matchings of 
facilities and sites with coefficients equal to the weights of the corresponding p.h.m.paths in (y,z). 

Theorem 35 The following holds true: Conv{QL) = Conv{Qj). 
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Proof. The theorem foUows directly from the combination of Theorem 
Corollary 



Theorem [271 and 



Corollary 36 The following mappings are bijective: 

1. Bi : ExtiQi) I — > 

2. B2 : Conv{QL) 1 — > ConviWi) 
Ext{WL) 

^ Conv{Wi) 
Ext{Conv{Wi)) 



3. S3 : Ext{QL) H 

4. Bi : Conv{QL) 

5. B5 : Ext{QL) ^ 



2.2 Model Objective 

Definition 37 (Objective function costs) Let (y,z) S Ql- V(i, j, n, u , A;, t) € F^, V(r, s) S : 
1 < r < s, the "cost" associated with Zui^virjkst is defined as: 

if r = 2; s = 3; i = v; k = j 
hu,i,t,s+i + /ii),2,t,s+i if r = 2; s > 4; i = u 
Oir + hirj^r+i + hirt,r+2 if 3<r<n-3;s = r + l;A; = j 



Cul,virjkst • — * 



Ojr + Oj,r+l + Ot,r+2 + hirj^r+1 + hirt,r+2 + /lj,r+l,t,r+2 

if r = n — 2; s = n — 1; k = j 
hirt,s+i if 3<r<n — 3;s>r + 2 
Otherwise 



Theorem 38 Let; 



l,virjkstZu,l,virikst- 



virjh 



(72) 



uGF vS-F iS-F rS-R: jSF fcGF sS-R: tS-F 



Then, i9(y, z) accurately accounts the cost of the perfect matching of the facilities and sites that is 
associated with (y, z) for all {y,z) £ Ext{Qi). 



(73) 



Proof. From Theorem | 

(y,z) G ExtiQi) =^ {y,z) G Qj 
Now, using Theorems [9l it can be verified directly that for (y, z) € Qi 

n—l n n 

'&{y,z)=Yl H hi^,r,is,s+ Yj "ir,r, where v E >I (?/, z) VreT. 

r=l s=r+l r=l 
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2.3 Overall Model 

Our overall linear programming model is as follows: 
Problem 39 {Problem LP) 

min {'d{y,z) : {y,z) G Ql} 

Theorem 40 The following statements are true of basic feasible solutions (BPS) of Problem LP 
and perfect matchings of the facilities and sites: 

1. Every BFS of Problem LP corresponds to a perfect matching of the facilities and sites; 

2. Every perfect matching of the facilities and sites corresponds to a BFS of Problem LP; 

3. The mapping of BFS's of Problem LP onto the set of perfect matchings of the facilities and 
sites is surjective. 

Proof. Statements (1401 1) and (j40[ 2) follow directly from Theorem I35( and the correspondence 
between BFS's of LP models and extreme points of polyhedra (see [H pp. 92-101]). Statement 
(j40[ 3) follows from the primal degeneracy of Problem LP (see |1H p. 32]). ■ 

Corollary 41 Problem LP and Problem QAP are equivalent; that is, Problem LP correctly solves 
the QAP. 

3 Numerical Implementation 

As indicated in section 1 and the abstract of this paper, we conducted a small numerical experi- 
mentation aimed at getting some rough idea about the computational performance of the model. 
In doing this, we implemented a streamlined version of the model (i.e.. Problem LP) in which con- 
straints (|12p-(jl4p were handled implicitly, and the upper bounds on the variables were omitted. We 
solved one set of five problems with 7 facilities, and one set of five problems with 8 facilities. Four 
(of the five) problems in each of the two sets were randomly generated, and one problem had all 
parameters (i.e., each inter-site distance, each inter-facility fiow, and each operating cost). equal to 
zero. In the randomly-generated problems, the inter-facility fiows and the inter-site distances were 
assumed to be uniform random numbers between and 300, and between 1 and 50, respectively. 
The facility operating costs were assumed to be zero in two of the randomly-generated problems 
in each set, and assumed to be random deviates on [0, 1000] for the remainder two problems in 
the set. Furthermore, one of the two problems with operating costs equal to zero had asymmetric 
inter-facilities fiows and asymmetric inter-sites distances, while the other two problems had sym- 
metric inter-facility fiows and symmetric inter-site distances. Similarly, one of the two problems 
with operating costs greater than zero had symmetric inter-facilities flows and symmetric inter-sites 
distances, while the other two problems had asymmetric inter-facility fiows and asymmetric inter- 
site distances. We solved each of the test problems using the simplex procedure on the primal and 
dual forms, respectively, and using the primal-dual interior-point ("barrier") method. The linear 
programming (LP) solver we used is the "Clp" set of routines of the COIN-OR open source library 

(S)- 
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The characteristics of the test problems and the computational results are summarized in Table 
1. The average computational time for the 7-sitc problems was 24.76 seconds, 59.44 seconds, and 
305.50 seconds of Sony VAIO notebook computer (1.83 GHz Intel Core 2 Duo processor) for the 
simplex procedure for primal form, the simplex procedure for the dual form, and the barrier method, 
respectively. The corresponding numbers for the 8-site problems, were 8,988.18 seconds, 3,064.86 
seconds, and 3,765.05 seconds, respectively. 

In general, it appeared that the computational performance of the model is not affected by the 
symmetry /asymmetry of the inter-site distances, nor by the symmetry/asymmetry of the inter- 
facilities flows, nor by the presence/absence of fixed facilties operating costs. Both the simplex 
method using the dual form, and the barrier method appeared to perform better than the simplex 
method using the primal form when the number of sites is increased. However, it also appeared that 
this could have been due to numerical stability problems. Overall, we believe the simplex procedure 
using the primal form may hold the greatest promise with respect to further developments aimed 
at solving larger-sized problems, based on the fact that it tended to examine a small number of the 
perfect matchings of the facilities and sites in general. 



Problem Details 


Simplex: Primal 


Simplex: Dual 


Barrier Method 






D(3) 


0(4) 


# of iter. 


CPU sec. 


# of iter. 


CPU sec. 


^ of iter. 


CPU sec. 


7 


A 


A 


Y 


19,167 


22.75 


16,434 


53.31 


17,048 


295.53 


7 


A 


A 


N 


24,610 


36.062 


21,189 


84.06 


22,326 


355.61 


7 


S 


S 


Y 


20,365 


27.14 


16,504 


50.84 


1,224 


272.91 


7 


S 


S 


N 


21.990 


34.52 


19.467 


71.20 


20,038 


441.92 


7 


N 


X 


N 


10.205 


3.31 


lO.GlO 


37.77 


5.284 


101.53 


Average 


19,274.4 


24.76 


16,840.8 


59.44 


13,184.0 


305.50 


8 


A 


A 


Y 


357,848 


10,837.66 


115,794 


2,763.28 


107,150 


3,589.19 


8 


A 


A 


N 


368,682 


11,023.23 


118,512 


3,036.61 


110,554 


3,760.24 


8 


S 


S 


Y 


379,912 


11,777.20 


132,518 


3,733.08 


124,685 


4,469.31 


8 


s 


s 


N 


:-!7().547 


11.2()(i.l!) 


131,022 


3,711.48 


124.897 


4.977.52 


8 


N 


N 


N 


32,831 


36.64 


67,559 


2,079.84 


69,498 


2,028.98 


Average 


303,164.0 


8,988.18 


113,201.0 


3,064.86 


107,356.8 


3,765.05 



1: Number of sites 

2: Flow type: "A"=asymmetric; " S" =symmetric; "N"=A11 flows are zero 

3: Distance type: " A"=asymmetric; "S"=symmetric; "N"=A11 distances are zero 
4. Operating Costs: "Y"=Not all equal to zero; "N"=A11 equal to zero 

Table 1: Summary of the computational experimentation 



4 Conclusions 

We have developed a linear programming model of the QAP. From a theoretical perspective, the 
proposed model provides a new affirmative resolution to the important issue of the equality of 
computational complexity classes P and NP. Although our proposed linear programming model has 
a polynomially-bounded size (O(n^) variables x O(n^) constraints), that size is still prohibitive 
with respect to practice. Hence, we believe developments aimed at solving larger-sized problems 
could be the subject of fruitful future research. 
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